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XMM-Newton observations of the supernova remnant IC 443: 
II. evidence of stellar ejecta in the inner regions. 
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ABSTRACT 

Aims. We investigate the spatial distribution of the physical and chemical properties of the hot X-ray emitting plasma of the supernova 
remnant IC 443, in order to get important constraints on its ionization stage, on the progenitor supernova explosion, on the age of the 
remnant, and its physical association with a close pulsar wind nebula. 

Methods. We present XMM-Newton images of IC 443, median photon energy map, silicon and sulfur equivalent width maps, and a 
spatially resolved spectral analysis of a set of homogeneous regions. 

Results. The hard X-ray thermal emission (1.4-5.0 keV) of IC 443 displays a centrally-peaked morphology, its brightness peaks 
being associated with hot (kT> 1 keV) X-ray emitting plasma. A ring-shaped structure, characterized by high values of equivalent 
widths and median photon energy, encloses the PWN. Its hard X-ray emission is spectrally characterized by a collisional ionization 
equilibrium model, and strong emission lines of Mg, Si, and S, requiring oversolar metal abundances. Dynamically, the location of 
the ejecta ring suggests an SNR age of ~4,000 yr. The presence of overionized plasma in the inner regions of IC 443, addressed in 
previous works, is much less evident in our observations. 

Key words. X-rays: ISM— ISM: supernova remnants — ISM: individual: IC 443 — pulsars: individual: CXOU J061705. 3+222127 



1. Introduction 

The Galactic SNR IC 443 (G 189.1+3 0) has an optical and ra- 
dio double shel l morphology dLeahvl I2004I: iReich .et all l2003t 
iBraun & Stroml 1 19861: iDuin & van der LaanUl975l) . Signatures 
of interaction with a dense and complex environm ent are visible 
at radio (e.g. ISnell et all 120051 : Ibenover 
lengths (Neufe ld & YUalll200alRho et al 



1977) and IR wave- 
2001, and references 



therein). In the northeast the shoc k front has been decelerated by 
the encounter with a H i cloud dRho et all 120011: iDickel et all 
1989; Denover, 1978). A giant molecular cloud, mapped by 
ICornett et all d 1977b . is located in the foreground, and is inter- 
acting with the remnant at several positions along the south 
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ern ri m (as schematically represented in Fig. 9 of Troiaet al 
2006). IC 443 is also a sou rce of y-rays (Esposito et al. 



1996; ISturner & Permed 1 1995b . Interestingly, in the same re- 
gion o f sky where OH (1720 MHz) maser emission is de- 
tected dHewitt et all|200dlHoffman et all 120031: IClaussen et all 
119971) . lAlbert et all d2007l) discovered a source of very high en- 
ergy y-ray emission, possibly originated by the interactions be- 
tween cosmic rays accelerated in IC 443 and the dense molecular 
cloud. 

Very interesting features and a peculiar morphology 
are displayed in th e X-ray band (Bvkov etall 120081 [2005; 
Bocchino & Bvkovi 120031 lAsaoka & AschenbacH 119941: 



iPetre et all fl988l) . A source of very hard X-ray emission, 
identified as a pulsar wind nebula (PWN), lies close to the 
SNR's southern rim and is possibly associa t ed with the remnant 
dGaensler et all 120061: iBocchino & Bvkovi 120011; lOlbert et all 
l200lh . 

On the basis of its X-ray properties, as observed by ROSAT, 
IC 443 has been classified as a mixed m orphology (MM) SNR. 
MM SNR main features were codified bv lRho & Petre ( 1991) as 
follows: 1) the X-ray emission is thermal and centrally peaked, 
with little or no evidence of a limb-brightened X-ray shell; 2) it 
arises primarily from swept-up interstellar material, not from 
ejecta; 3) a flat temperature radial profile, and a constant or in- 
creasing density toward the remnant center; 4) a physical asso- 
ciation with dense molecular clouds. 

XMM- Newton high resol ution observations of IC 443, pre- 
sented by iTroia et al] d2006l) . resolved a limb-brightened shell 
structure in the very soft X-ray band (0.3-0.5 keV), and reported 
a first evidence of Mg, Si and S enhanced abundances. In addi- 
tion, XMM-Newton and Chandra observations of several other 
MM SNRs unveiled very steep X-ray brightness profiles and 
multi-c omponent X-ray-emitting plasma with enhan ced abun- 



Wan g et all 119921: IPetre et all Il988l) . The bulk of its X-ray 



emission is thermal, fairly described with a two components 
ionization equilibrium model of temperatures ~0.2-0.3 keV 
(hereafter "c old component") a nd >1.0 keV (hereafter "hot 
component' '; ITroia et a l. 2006; lAsaoka & Aschenbachl 119941: 



dances dLazendic & Slanel 120061 IShelton et alll2004l) . More re- 
cent results suggest indeed a non negligible presence of metal 
enrichment gradients and stellar ejecta inside at least ~50% of 
all known MM SNRs. 

Strong Si an d S lines in the IC 443 X-ray spectrum were 
fir st reported by IPetre et all d!988l) . and then studied in detail 
bv lKawasaki et al.1 d2002l) . In the latter work the X-ray emission 
from IC 443 was described by a two temperatures model, with 
a 0.2 keV component in ionization non-equilibrium (NEI) and a 
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1.0 keV component in collisional ionization equilibrium (CIE). 
Narrow Gaussians were used to reproduce the most prominent 
lines. The elements abundances were all below the solar value, 
as expected for X-ray emission from the shocked ISM. Fro m the 
relative strengths of Si and S lines, iKawasaki et al.l d2002l) de- 
rived an ionization temperature of -1.5 keV, significantly above 
the 1 .0 keV continuum temperature. Those results led to the hy- 
pothesis of an overionized thermal plasma. 

Using XMM-Newton observations of the SNR IC 443, we 
have already presented a detailed analysis of the soft X-ray ther- 
mal emission, mainly o riginated from the sho cked ISM (see 
iTroia et al.ll2006l) . Unlike Kawasaki et al. I (120021) our spectral re- 
sults gave a soft X-ray component (kT~0.2-0.3 keV) near the 
equilibrium condition (r~10 12 crrT 3 s) in most of the analyzed 
regions. The hot component (kT>1.0 keV) was in full equilib- 
rium and showed a high metallicity. 

The present work focuses on the properties of the hot X-ray 
thermal component, and is aimed at 1) addressing the high metal 
abundance plasma inside the remnant and its ionization stage; 2) 
checking whether the abundance pattern of the hot X-ray emit- 
ting plasma is consistent with a core-collapse SN as progenitor, 
as already hinted by the presence of the PWN. 

The paper is organized as follows: observations are described 
in §12 results from imaging and spectral analysis are presented 
in §[3] In §|4]we discuss their implications for the remnant ion- 
ization stage ( §14. U . the SN progenitor ( §I4.2| |. and the age of 
IC 443 ( jj4.31 l. A summary of our conclusions is presented in §0 
Throughout the paper, the quoted unc ertainties are at the 90% 
confidence level (Lamp ton et aU [l976). unless otherwise stated. 

2. Observations and Data Reduction 

We analyzed four pointings XMM-Newton observations of 
IC 443, performed during the Cal/PV phase on 2000 Septe mber 
25-28 and previously d escribed in Troi a et aTl (2006) and 
iBocchino & Bvkov (2003;). We added a more recent and very 
deep XMM-Newton observation, taken on 2006 March 30 
(ObsID 0301960101, P.I. F. Bocchino), and centered on the 
southeastern side of the remnant at R.A.(J2q00)-0 6 h 1 8 m 04 s , 
Dec(J2000)=+22°27'33" (see also lBvkov etafll2008l) . 
Our work is focused on the analysis of EPIC MOS data. We used 
the Science Analysis System software (SAS , version 7.0) for data 
processing and xspec v. 1 1 .3 ( Arnaud, 1996) for spectral analysis. 

Data were screened to remove spurious events and time 
intervals with heavy proton flaring. After the particle back- 
ground removal, the total exposures for the MOS cameras 
(MOS1+MOS2) were -170 ks and -140 ks for the 2000 and 
2006 data sets, respectively. 

In order to estimate the Galactic background we used a set of 
blank sky observations performed as part of the Galactic Plane 
Survey (GPS, P.I. Parmar). We verified that, in the energy band 
(0.3-5.0 keV) selected for the analysis, our results are not sensi- 
tive to the background subtraction (either GPS or hig h Galactic 
latitud e blank fields). We refer to our previous paper dTroia et al.L 
2006) for more details about the observations and the data reduc- 
tion techniques we used. 

3. Results 

3. 1 . X-ray Images 

Fig. Q] shows the composite image of IC 443 in the hard 1.4- 
5.0 keV energy band, where the contribution of the cold X-ray 
emitting plasma is negligible (see lTroja et al . 2006). The image 
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Fig. 1. X-ray count rate image of IC 443 in the 1.4-5.0 
keV energy band. Units are ctss^bin -1 . The image was 
adaptively smoothed using <x m ,„=5" and cr max -20"; the bin 
size is 5". Contour levels correspond to 0.4, 1.0, 2.7, and 
5.8xl0~ 2 ctss" 1 arcmin -2 . Black crosses mark the position of 
the 12 hard X-ray point sources identified bv lBocchino & Bvkovl 
(2003). The dashed line traces the symmetry axis of the X-ray 
emission in the hard band; the black arrow indi cates the puta- 
tive di rection of motion of the PWN, according to lGaensler et al.l 
J2006h . 



is background subtracted, exposure and vignetting corrected. We 
re ported the position of the 12 hard X-ray point sources detected 
bv lBocchino & Bvkovl (120031) . 

The hard band morphology does not resemble IC 443 emis- 
sion in other bands (optical, radio, IR and soft X-ray) and, on the 
plane of the image, it appears inner concetrated with respect to 
the cold X-ray emitting plasma. 

The hard X-ray thermal emission is dominated by an elon- 
gated axis-symmetric structure, having its surface brightness 
peak in the northeastern (NE) quadrant. The symmetry axis, 
traced by the dashed line in Fig. [T] does not match with the ple- 
ri on wind nebula (PWN ) direction of motion (-130° N through 
E. lGaensler et alj|2006l) . as indicated by the arrow in Fig.Q] 

It is worth noting that the X-ray surface brightness is higher 
(2 > 0.012 cts s 1 arcmi n -2 ) in the subshell A (following the 
nami ng convention of iBraun & Stroml 119861; see also Fig. 5 
of ITroia et alj [2006), which is associated with the presence 
of the neutral cloud in the NE, while it appears strongly re- 
duced (2 < 0.012 cts s -1 arcmin~ 2 ) in the western part (sub- 
shell B), where the blast wave is expanding in a more rarefied 
(no~0.25 cirT 3 ) and homogeneous medium. 

3.2. Median Energy Map 

XMM-Newton great sensitivity and collecting area enabled us to 
create a high resolution image of the median photon energy £50 
in the 1 .4-5.0 keV band. In the selected energy range the median 
energy is mostly influenced by the parameters of the hot X-ray 
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component, being unaffected by spectral variations of the colder 
one. For a plasma in ionization equilibrium and with solar abun- 
dances (model MEKAL in xspec) we verified that the £50 dis- 
tribution represents a powerful tracer of temperature variations, 
especially forkT>l keV. 

Fig. [2] shows the median energy map of the X-ray emitting 
plasma in the 1.4-5.0 keV band. For a comparison, we overlaid 
the brightest hard X-ray emission contours of Fig.[TJ The white 
dashed line guides the eye to identify a ring-shaped structure, 
which will be discussed in j 



Table 1. Energy bands selected for the equivalent width images. 



The median photon energy image clearly shows a north-to- 
south elongated structure, spatially corresponding to the bright- 
est emission features in the hard energy band (contours in Fig. [2]). 
The region of high median energy (£5o>1.9 keV) along the 
southeast edge can be partially attribut ed to the emission of th e 
PWN and of the other point sources dBocchino & Bykov, 2003|), 
but it appears to have a much larger extent, suggestive of an ex- 
tended hard emission. The nature of the hard X-ray source 1 SAX 
J06 1 8.0+2 227 (also marked in Fig.lp has been extensively in- 
vestigated dBvkov et all 12008, 2005; B occhino & Bvkovi feOOO) 
and interpreted in the framework of interaction of a molecular 
cloud with a fast moving SN ejecta fragment. The study of the 
diffuse hard emission revealed in Fig. [2] is beyond the scope of 
this work, however we note that it may have a similar physical 
origin, probably arising from the interaction of the SNR with th e 
structured cloud (van Dishoe c¥et~aTlll993l:lB"urton et all [l990). 
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3.3. Equivalent Width Images 

IC 443 X-ray emission is primarily thermal and its spectrum is 
dominated by prominent atomic emission lines. Its global spec- 
trum can be described with a phenomenological model, con- 
sisting of two bremsstrahlung components for the continuum 
(kT~0.2 ke V and kT~ 1 . 1 ke V) plus Gaussian features reproduc- 
ing the line emission. 

We explored the spatial distribution of the hard X-ray lines 
emission by constructing equivalent width (EW) images, often 
used as powe r ful tracers of stellar ejecta fragments i n SNRs 



dHwang et all 120001: 



Cassam-Chenai et al 



Gotthelf etall [200Tt IParketaTl [2003a; 



20041 iMiceli etafl I2006T) . 



We selected an energy band for each emission line, com- 
bining the Hea and Lya emission of each element (e.g. Si XIII 
and Si XIV) in a single image to improve the statistic. We 
chose the energy band that optimizes the signal to noise 
S/N) ratio according to ou r best-fit phenomenological model 
ICassam-Chenai et alll2004l) . 

X-ray emission in the line energy band needs to be cor- 
rected for the underlying continuum emission. The correction 
is made pixel by pixel because the continuum may be spatially 
distributed in a different way from the line emission. In order 
to account for the continuum contamination, we selected two 
continuum bands located on either side of the line energy band. 
Their width was chosen to optimize the S/N ratio, according to 
the recipe of ICassam-Chenai' et al.l (120041) . 

This method was successfully applied to silicon and sulfur 
lines. We excluded the H-like Mg line at 1.47 keV since the un- 
derlying continuum emission cannot be correctly estimated. In 
this particular case, the low energy continuum band (~1.1 keV) 
is dominated by the soft spectral component and still affected 
by the strong absorbing column gradients, due to the foreground 
molecular cloud. 

Line and continuum band passes, selected for each spectral 
line of interest, are listed in Tab. [TJ 

Once an image in each energy band was generated, we 
estimated the continuum e mission in the line band (see e.g. 
Cassam-Chenai et al. 2004). In each pixel we computed the un- 
derlying continuum as a linear combination of the source emis- 
sion in adjacent continuum bands: 



I c =Aa:i: + (1 -X)c£ t c 



(1) 



Fig. 2. Median photon energy map in the hard band (1.4- 
5.0 keV). Units are eV. Hard X-ray emission contours, cor- 
responding to 2.7 and 5.8xl0~ 2 cts s" 1 arcmkr 2 , are overlaid. 
Noisy (S/N<3) regions inside the FOV, delimited by white solid 
line, were masked. The bin size is 20"and the smoothing width 
is cr=20". The two dashed semi-circles guide the eye to notice 
a ring-shaped structure, which is more evident in the EW maps 
shown in Fig. 3. The location of some sources, partially respon- 
sible for the high median energy in the SE region, is indicated 
with green arrows. 



where I c is the resultant continuum image, and /* are the 
background subtracted images in the high and low contin- 
uum bands, respectively. The normalization coefficients, a*, 
were determined by the ratio of the model predicted counts 
and the source measured co unts in the continuum bands, as in 
ICassam-Ch enai et al. (2004). The weighting factor A was tuned 
so to minimize statistical fluctuations. 

Narrow band images were created with a bin size of 30" to 
provide at least 10 counts even in the fainter SNR regions. 

In order to generate EW images, the image in the line band 
was background and continuum subtracted and then divided by 
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Fig. 3. Equivalent width images of Si (fo/7 panel) and S lines 
(bottom panel). Units are eV. Noisy regions were masked, as de- 
scribed in the text. The thin white line mark the FOV. Contour 
levels corresponds to 450, 550, 650 eV in the Si EW image. We 
also overlaid the three regions (named IN, RING, and OUT) se- 
lected for spectral analysis. The dashed boxes overlaid on the S 
EW map show the " North" and "Center" regions, the same of 
iKawasaki et ail d2002l) . and discussed on jH3.4. II 



the estimated continuum image I c . Since this operation ampli- 
fied statistical fluctuations, we previously applied an adaptive 




Energy [keV] 

Fig. 4. A comparison of the three X-ray spectra at high energies, 
clearly showing brightest lines in the "ring" region respect with 
the external ("out") and the internal ("in") regions. 



smoothing with 30" < <x <120"and a signal to noise ratio equal 
to 5. We finally set the EW to zero in those pixels with a very 
low estimated conti nuum (<5% than the maximum value, see 
e.g. lParketal1l2003bl) . 

Si and S EW images are presented in Fig. [3] revealing an in- 
triguing circular feature, 10'in radius, in the southern part of the 
SNR. For a morphological comparison we superimposed Si EW 
contours on the images, corresponding to 450, 550 and 650 eV. 
The two images look very similar. The contours of high Si EW 
match well with the bright S EW features. In particular, we did 
not find any strong evidence of elemental stratification. Both the 
EW maps show similar features to the median energy map: an 
elongated bright structure in the NE quadrant and a circular ring 
around the PWN region. A rather similar arc-like feature has 
been o bserved in the Si EW map of the SNR N49 dPark et all 
l2003al) and interpreted as the boundary of hot plasma, re-heated 
by the reverse shock. The Si arc in N49 is not around the central 
compact object but, as in the case of IC 443, it is located near the 
site of interaction with a molecular cloud complex, strengthen- 
ing the hypothesis of an intense reverse shock, propagating back 
into the remnant. 

The great advantage of EW images is that they are not signif- 
icantly affected by emission measure variations. However they 
do also depend on the plasma temperature and ionization age. In 
order to disentangle multiple effects, we proceeded with direct 
spatially resolved spectral analysis, using the median photon en- 
ergy and the EW images as guiding criteria of region selection. 

3.4. Spectral Analysis 

Our aim was to confirm the bright X-ray features visible in Fig. [2] 
and Fig. [3] and give a quantitative characterization of their spec- 
tra. We extracted three EPIC MOS spectra, selected as represen- 
tative on the basis of the median energy map and the EW images. 
We selected homogeneous regions, marked in Fig. [3] character- 
ized by small fluctuations (<<5%) in the median photon energy 
and in the line EWs. Two of them, named "in" and "out", have 
low median energies, low EWs values, and lie inside and outside 
the ring structure, respectively. The third region, named "ring" 
and belonging to the ring structure, is particularly bright both in 
the median energy and in the EW images. 
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Table 2. Spectral parameters of the three analysed regions. 
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Fig. 5. EPIC MOS spectra of IC 443, selected as representative 
of the hot X-ray plasma properties, and their best-fit models 
overlaid. Panels a, b, and c correspond to regions "in", "ring", 
and "out" in Fig. [3] respectively. The lower panel in each plot 
shows the residuals from the best-fit model. 



Source and background events list were corrected for vi- 
gnetting effects with the task evigweight of the SAS package, 
and then spectra were extracted. We used the relevant on-axis 
ancillary response file and response matrices generated with the 
SAS task rmfgen. Spectra were grouped to a minimum of 25 
counts per bin and the x 2 statistics was used. Fig. [4] show the 
three spectra at high energies, where striking differences in the 
properties of the hot X-ray emitting plasma are clearly visible. 
Prominent Mg, Si, and S lines characterize the "ring" spectrum, 
whilst line emission is much fainter in the other two regions. 
Fig.|4]nicely confirms that the identified features correspond to 
different physical conditions of the emitting plasma. 

We modeled the spectra over the 0.5-5.0 keV en- 
ergy band with a NEI plus a CIE t hermal componen t 
(VPSHOCK+VMEKAL in xspec), as in [Troia et alj (I2006h . 
MOS1 and MOS2 spectra were simultaneously fitted, leaving 
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Ej (keV) TTO" 

SiEW(eV) 310 

S EW (eV) 280 

Cold component 

N H (10 22 cm- 2 ) [061] [538] 

kT, (keV) 0- 44 !ora 0.43±0.02 

T(10 12 cm- 3 s) l-7!oo6 1-8+0.05 

O/O 0.78 :;:" 0.68 ;; ; ; 

Ne/Ne G O.'H ;; ! :: 0.81!°f 

Mg/M go 0.451°;* 0.35±0.05 

Si Si. 0.26 0.50!™ 

Fe/Fe Q 0.10 0.20±0.03 

EM,,(10 56 cm" 3 ) 4.1 

Hot component 

kT, (keV) 05^ 

Mg/Mgo 2.6!^ 

Si/Si 1.1+0.3 

S/S G 0.7±0.5 

Fe/Fe Q <0.13 

EM ; ,(10 56 cm" 3 ) 1.1 

^ 2 /dof 420/324 
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1.65 
290 
130 



5.7 



[0.50] 

32 +003 

"• J -0.02 

2 5 +0 ' 5 
0.30+0.10 
35 +01 ° 
0.35±0.05 

<0.15 
0.22±0.05 
8.7 



1 4C+O.U8 

5 +3 

J - -0.9 

2.4+0.7 

9 q+0.6 
^ -0.7 

<0.03 
0.4 
497/343 



0.86+0.05 
[1] 

67 +019 
22 +012 
<0.23 
0.4 
412/302 



the normalization between the two instruments as a free param- 
eter. We included a 5% systematic error term accounting for the 
uncertainties in the calibratiorfj The absorbing column Nh was 
estim ated by using the median photon energy in the soft band 
(cf. lTroia et alj 2006, Fig. 6), and held fixed at the derived value. 

Spectra and best-fit models are shown in Fig. [5] Table ^re- 
ports the spectral properties (median photon energy E50, Si and 
S EWs), and the best-fit spectral parameters for each region. 



3.4.1. Overionization 

Taking advantage of the large XMM-Newton effective area at 
high energies, we checked the overionization claim made by 
Kawasaki et al. (2002, 2005). For this purpose, we extracted two 
spectra from the "No rth" and "Center" regions (see Fig. 2 of 
iKawasaki et al.ll2002l) . and performed the same comparison be- 
tween the ionization and the electron continuum temperatures. 

We selected the 2.2-5.0 keV energy range and we mod- 
eled the spectrum with a thermal component (model VMEKAL 
in xspec), fixing the S abundance to zero. Three narrow gaus- 
sians (centered at 2.46, 2.62 and 2.88 keV) were used to model 
the sulfur line emission. The intensity line ratio of the H-like 
and He-like ions (S xvi/S xv) is a function of the S ionization 
temperature, derive d using the APED line emissivity database 
(ISmith et al.L[200lb . 

The absorption column densities were held fixed to the 
values of 0.60xl0 22 cirr 2 for the "North" spectrum and 
0.68xl0 22 cm 2 for the "Center" sp ectrum, as derived from the 
Nh/Eso calibration plot (Fig. 10 in iTroia et af] |2"006). We also 



1 The current EPIC MOS status of calibration is reported at 
http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0018.pdf 
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Fig. 6. Confidence level ^contours (68%, 90%, 99%) of the line 
ratio (S xvi/S xv) and electron temperature (kT e ) for the "North" 
and the "Center" region. On the right Y-axis is reported the ion- 
ization temperature kT z . The dashed line indicates the CIE state. 

tested the absorption c olumn value of 0.74xl0 22 cm 2 , used by 
iKawasaki et al.l (|2002|) . but it does not affect the final result as 
expected at these high energies. 

Fig.[6]shows the 68%, 90%, and 99% x 2 contours in the pa- 
rameters space S xvi/S xv- kT e . The CIE is represented by the 
dashed line, tracing the condition kT e =kT z . In the "North" re- 
gion we found evidence of an underionized plas ma, in stark con- 
trast with the findings of Kawasa ki et al.l d2002l) ; in the "Center" 
region the derived ionization temperature is slightly higher than 
the the electron temperature,fhough consistent with it at the 90% 
confidence level. Therefore, a detailed analysis of the ionization 
state does not directly support the presence of an overionized 
plasma in the "North" region, and shows only a marginal (1 cr) 
evidence of overionization in the innermost "Center" region. 

4. Discussion 

4.1. Presence of hot metal-rich plasma in IC 443 

Spatially-resolved spectral analysis confirms that the structures, 
revealed by the median photon energy and the EW maps (see 
33.2l and 33.31 respectively), are associated with a high tempera- 
ture, high metallicity X-ray emitting plasma. 

Our results, listed in Tab. [2j show that the cold plasma is 
characterized by a temperature of 0.3-0.4 keV, a high ionization 
timescale (>10 12 cirT 3 s), and solar or undersolar abundances, as 
expected if the emission arises from the shocked ISM. The hot 
plasma is fully equilibrated. It has lower metal abundances in the 
"in" and "out" regions, in good agreement with previous finding s 
for the hot component (IKawasaki et all 120021 : iPetre et all [1988). 
The temperature T;, and the abundances of Mg, Si, and S in the 
"ring" region are much h i gher, m ore similar to regions 2, 6, and 
9 analyzed in lTroja et alj d2006l) . As shown in Fig. [2] and Fig. [3] 
in the NE the hottest plasma appears concentrated in the inner- 
most regions, whilst in the south it display a ring-shaped spatial 
distribution. Our spectral results fairly agrees with this picture, 
detecting a high temperature and a substantial metal-enrichment 
in correspondence to the "ring" region. The measured overabun- 
dant Mg, Si and S are suggestive of SN ejecta emission. 

IC 443 is a middle-aged remnant, thus the presence of ejecta 
at its evolutive stage is not expected. Models predict that a great 
part of the ejecta thermalized with the surrounding medium, 



being no longer visible. However, XMM-Newton and Chandra 
high-resolution observations show a growing number of middle- 
aged remnants (~10 4 yr) where metal rich ejecta are still de- 
tectable an d significantly contr ibute to t he X-ray emission, as in 
DEML71 dHughes et all 120031) . N49B dPark et all l2003al) . and 
0103-72.6 dPark et a l.. 2003b). Moreover, there are several cases 
among MM SNRs in which a soft component, associated with 
shocked ISM emission, and a hot metal enriched compo nent 
have been resolved. For instance, Mavromat akis et all d2004l) de- 
scribed W63 thermal emission with a two temperatures plasma 
(kTi~0.2 keV and kT2~0.6 keV), showing overabundant Mg, Si 
and Fe i n a bright bar-like region near the geometric center of 
the SNR. lLazendic & Sland (2006) resolved a two temperatures 
thermal plasma in the inner regions of CTB 1 , a prototypical MM 
SNR, and they showed that the hard component (kT~0.8 keV) 
requires an oversolar Mg abundance. 

IKawasaki et al.l (I2002L [2005) addressed the issue of overi- 
onization in the hot plasma of MM SNRs, finding that in 2 
objects (namely IC 443 and W49B), out of a sample of 6, 
some ions are indeed overionized. Such ionization state of the 
plasma has many important implications, since it would indicate 
rapidly cooling plasma due, for instance, to thermal conduction. 
A thermally conductive model has been invo ked to explain the 
centrally peaked mo rphology of MM SNRs dCox et all Il999fc 
IShelton et all Il999l) . and overionization would provide a nice 
consistency check. 

We investigated whether overionization conditions are con- 
sistent with our XMM-Newton data, performing our check i n the 
same regions of IC443 analyzed by IKawasaki et ail d2002l) . In 
the "North" region, we obtained that overionization is not con- 
sistent with the results for the S ions. In the "Center" region, 
overionization is still possible, though not required at the 90% 
confidence level, for the same ions. We point out that this region 
is characterized by a wide spread of EWs in both the Si and the 
S maps (AEW/<EW>~40%), therefore the analyzed spectrum 
originates from chemically inhomogeneous regions. A check of 
the ioniz ation state would re quire more uniform conditions. For 
instance, iMiceli et a i] d2006l) found that in the central core of the 
SNR W49B there is no compelling evidence of overionization. 
Unfortunately, our attempt to analyze smaller but more homo- 
geneous regions was affected by the low statistics, and did not 
significatively constrain the spectral parameters. 



4.2. Constraints on the SN progenitor 



We presented in § l3.2H3.3l high resolution maps of some proper- 
ties of the hot plasma (median photon energy, and EWs of Si and 
S lines). We further investigated and quantified the properties of 
the X-ray emitting plasma through a spatially-resolved spectral 
analysis (§ 13.41 ). Those results allowed us to infer the spatial dis- 
tribution of the stellar ejecta, and to estimate their abundances 
pattern and masses in the ring structure, which are strong clues 
to the nature of the SN progenitor (e.g. [Rakowski et al., 2006; 
IMiceli et ail 120061: [Park et alll2004l) . and the age of the remnant 
(see § 1431 

Tab. [3] reports the best fit metal abundances of the "ring" re- 
gion relative to Si (col. 1), the abundances ratios expected for 
core-collapse SNe with different progenitor masses (col. 2-4; 
IWooslev & Weaver! 1 19951) . and for Type la SNe according to 
different deflagratio n models (col. 5-6; [Rakowski et al. 2006; 
iBadenes et al.l l2003). Bearing in mind the uncertainty in the pre- 
dictions of nucleosynthesis yields, our comparison is aimed at 
discriminating between Type la and Type II SN explosions. Any 
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further inference, for example on the mass of the progenitor, 
would be highly dependent on the adopted model. 

The very low Fe abundance can not meaningfully constrain 
the progenitor explosion, and it can be explained if the innermost 
ejecta layers have not been shocked yet. The high Mg abundance 
allowed us to disfavor a Type la SN, which primarily produces 
Fe-group elements and a negligible fraction of lighter metals. A 
Type II SN has been previously suggested as the remnant pro- 
genitor, only on the basis of its location in a star-forming region 
and the putative association with a PW N (Gaens ler et all 12006: 
iBocchino & Bvkovi 120011: lOlbert et all 120011) . Our results fur- 
ther support and strengthen a core-collapse origin for IC 443. 

In order to estimate the total mass of the ejecta, we as- 
sumed that the ejecta are confined within a shell of radius 
R e j~4.5 pc, and width Ar~l pc, for a total emission volume 
V s / lc //~7.7 x 10 57 cm 3 . Hence the hot plasma fills the entire vol- 
ume V S heii (filling factor //, = 1). The above assumption might not 
be valid for the Mg ejecta, since it is based on the EW maps of 
Si and S line complexes. Such high-Z elements are thought to be 
synthesized by explosive O- or Si-burning in the core of a mas- 
sive star, while lighter metals mostly reside in the outer layers. If 
the stellar stratification persists after the explosion, core-collapse 
nucleosynthesis models predict the Mg ejec ta to be preferen- 
tially located outside the Si shell (cf. Fig. 1 of Thielemann et al., 
1996), thus occupying a larger volume. The possibility of mix- 
ing or overturning of the ejecta layers (Willin gale et all [20021: 
iHwang et all [2000) add further uncertainty. Therefore the esti- 
mated mass of Mg should be taken with care. 

For the case where the electrons are primarily from H/He 
(« e =1.2«#), we derived a hot component gas density of 
«//=(0.62+0.10) cirT 3 in our "ring" region, which intercepts 
a line of sight of ~2 pc. By assuming the same line of sight 
for the inner and outer regions, we estimated a density of 
wh=(0.60+0.10) cm -3 , and«H=(0.96+0.16) cm -3 , respectively. 

We extrapolated the density derived in the "ring" region 
to the whole shell of ejecta, thus estimating the total mass of 
hot metal-rich plasma, M tot ~4M . We calculated a lower limit 
on t he mass of ejecta M z as M z /M tot > A z x (M z /M tot ) 
(e.g. lLazendic et all 120051) . where M tot is the total mass of 
the hot X-ray-emitting plasma, A z is the measured abundance 
for the element Z, and (M z /M tot ) is the standard solar mass 
fracti on (6.4 X 10' 4 for Mg, 7 Ox 10 4 for Si, and 3.6x10 4 
for S: lAnders & Grevessdfl989l) . We obtained: M Mg >0.013M G , 
M Si ^0.007M o , M s >0.003M o . 

An upper limit on M z was derived by the assumption of 
"pure" Z-ejecta, i.e. stellar fragments in which all the elec- 
trons are provided by the ionized specie Z, n f ~K z n z . This 
is a crude approximation of the ejecta mass, but not unreal- 
istic as this kind of objects have already been observed in 
other SNRs, as the Si-rich knots in Cas A (Laming & Hwang, 
120031) . We considered the contribution from both the ionization 
states of Si and S (K si =12.5, K s = 14.5), but only the highest 



Table 3. Relative abundances in the ring region: [X/Si]/[X/Si] 



Ratio 


IC443 


2OM fl 


25 M Q n 


30M o fl 


PDDe" 


DDTe 6 


Mg/Si.. 


2 1 +L4 

z..i_ 07 


0.16 


0.45 


1.79 


0.0017 


0.025 


S/Si 


1.0+0.4 


1.28 


0.96 


0.24 


1.5 


1.4 


Fe/Si... 


<0.018 


0.88 


0.13 


0.15 


0.89 


0.91 



a IWooslev & Weaver! jl995l) . * Pulsating d elayed detonation (P DDe) 
and delayed detonation (DDTe) models from Badenes et al. (2003). 



2.0 



1.5 



0.5 



0.0 



1 

s=0 uniform medium 


hlastwave shock 


n=0 uniform ejecta 




7 Subshell B Radius 




Ring Radius — ~~ 


.-4000 yr 

/ x. reverse shock 





Fig. 7. Shock positions for blast-wave (thin solid line) and re- 
verse (thick solid line) shocks in r emnants of uniform ejecta and 
evolving in an uniform medium dTruelove & McKeeL Il999h . A 
comparison with the observed values for IC 443 (ri~ll pc for 
the subshell B, and r2~4.5 pc for the ring of hot gas, tracing the 
forward and reverse shock location, respectively) gives a SNR 
age of 4 x 10 3 yr. 



ionization state of Mg (KM g =ll), as results from the best fit 
model (see Fig. [5] middle panel). We obtained: MM g ^O.55M , 
M si <O.3OM o ,M s <O.15M . 

It is worth noting that according to IWooslev & Weaver! 

(1995) the predicted masses of Mg, Si and S ejecta lie in the 
ranges between 0.01-0.4 M , 0.02-0.4 M , and O.O1-O.2M , re- 
spectively. These values well agree with the ranges derived from 
our analysis, which can be considered meaningful constraints on 
the ejecta masses in the hot shell. 

4.3. Age of IC 443 

We examined the case of a SNR expanding in an uniform 
ambient medium (s=0) and with an uniform ejecta distribu- 
tion (n=0). Similar results are found when considering a power 
law ejecta distribution (n=2). We took our estimate of the pre- 
shock density in the eastern regions (subshell B), no~0.25 cm -3 , 
and the mass of the hot plasma calculated in 34.21 For typi- 
cal SN energies (~10 51 erg ) , the l ength scale and the time scale 
of the ITruelove & McKeel d 19991) solution are R ch ~8 pc and 
t c h ~2.3 x 10 3 yrs, respectively. 

Fig. [7] show the predictions of ITruelove & McKeel (Q999) 
compared with the observed parameters of IC 443: a forward 
shock radius of 1 1 pc (shell B radius) and a reverse shock ra- 
dius of 4.5 pc (ring radius), correspond to an evolutive age of 
~400 yr, much youn ger than some previous estimates. In partic- 
ular, [Chevaliea Cl999|) derived an age of 3 x 10 4 yr by assuming 
that the remnant has been evolving in a very dense environment 
(~15 cirT 3 ) since its early sta ges. Our result is instead more sim- 
ilar to the age estimated by iPetre et al] dl988l) . who considered 
an exponential density distribution for the circumstellar environ- 
ment. The analysis of the soft X-ray emission fr om IC 443, and 
the derived density profile (cf. lTroia et all 12006, Fig. 7) suggest 
that the blast wave propagated in a much lower density medium, 
and only recently has encountered a sharp density gradient in 
the NE (subshell A), impacting onto a neutral hydrogen cloud 
dRho et al.Ll200ll:lDenoveii [1978I) . 
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As derived in 34.21 the remnant likely originated from the 
collapse of a massive star. Hence, during the progenitor life- 
time strong stellar winds could have blown away the surround- 
ing molecular clouds, formi ng a low densi t y cavi ty, where the 
remnant is now evolving. Brau n & Stroml dl986l) first argued 
that the multiple shell morphology of IC 443 might be the re- 
sult of such pre-SN interaction, and that the blast wave is im- 
pacting on the pre-existing wind-bubble walls. A detailed kine- 
matical study of the o ptical filaments in IC 443, carried out by 
Meaburn et all d!990h. gave furt her support to the scenario pro- 
posed by iBraun & Stroml d 19861) . The authors detected high ve- 
locity features along their slit Pos. 9, at the edge between sub- 
shell A and B, which can be easily explained if the emitting ma- 
terial, accelerated by the passage of the shock, is seen edge-on. 
The resultant geometry of the optical filaments is fully consis- 
tent with a n SNR evolving with in pre-existing interconnected 
cavities (cf. lMeaburn et all 19901 Fig. 5). 

Since the medium surroundin g massive stars is considerably 
modi fied by pre-SN winds (e.g. IChu et al.L 120031 : iDwarkadas. 
120051) . we verified whether our estimate of the age of IC443 
significantly depends on the assumption of a uniform ambient 
medium. For this purpose, we modeled the evolution of a SNR in 
an isothermal medium with a power-law density profile pocr~ 2 , 
by performing a numerical hydrodynamic simulation. We used 
the FLASH code (Fryxell 2000) and we set up a spherically 
symmetric explosion of 3 M of ejecta with a total energy of 
10 51 erg. The initial velocity profile of the ejecta increases lin- 
early with their distance from the center, and their density profile 
decre ases exponentially (in agreement with IWang & Chevalier! 
120021) . In this case p s =M w /(4w w )~1.76 x 10 14 gcm _1 , where 
M w is the mass-loss rate of the wind and v w is its speed, and 
both are assumed to be constant. The location of the forward 
and reverse shock fronts obtained in the simulation suggests a 
slightly younger age (~3000 yr). Our basic assumption, i. e. that 
the radial separation between the forward and the reverse shock 
is the result of the SNR evolutive stage, does not hold if the for- 
ward shock in the NW has been strongly decelerated in the cavity 
walls. Modifications of the shock velocity, due to the interaction 
of the blast wave with the confining walls, led to systematically 
underestimate the SNR age. In that case, the inferred value of 
4000 yr only provides a lower limit. However, in such a scenario, 
a close correlation between the two radii, as shown in Fig. [7] is 
not expected, as the forward shock radius would mainly depend 
on the properties of the cavity wall. We therefore conclude that 
the morphology of the X-ray emitting plasma in IC 443 robustly 
suggests an age <10 4 yr. 



4.3.1. Implications for the PWN 

A detailed study of the morphological and spectral properties 
of the PWN, and its central neutron star (NS; J0617+2221) has 
been presented bv lGaensler et al.l d2006l) . They inferred a black- 
body temperature kT°°~100 eV for the NS, and argued that the 
properti es of the thermal radiation from the NS are consistent 
with the lChevarieridl999l) age of 30,000 yrs. B y using the theo- 
retical cooling curve of lKaminker et alj d2002l) . we checked that 
our new estim ate of the SNR age and the temperatures range 
measured by Gaensleret al. (2006) are still consistent with the 
cooling of a 1.35-1.45 M NS with a lp proton superfluid core. 
Thus, the present knowledge of NS cooling curves does not 
allow us to put tighter constraints on the NS age, that would 
strengthen the associati on between the rem nant and the PWN. 
The NS temperature of Gaensler et al. (2006) and a putative age 
<10 4 yr are in agreement with observational results from other 



young NSs: PSR Jl 119-6127 (-1.6 kyr) RX J0822-43 (-3 .7 
kyr) IE 1207-52, and PSR J 1357-6429 (~7 kvr: IZavlinll2007bO: 
iKaminker et al.ll2002l) . 

The location of the ring of hot metal-rich plasma in our 
Fig. 2 and Fig. 4, just around the PWN, is an additional argu- 
ment in support of the physical connection between the PWN 
and IC 443. On the plane of the image (see e.g. Fig. 2) only the 
south side of the PWN seems to be in contact with the reverse 
shock, suggesting that the PWN may be in the reverse shock or 
subsonic expansion stage, bu t not in the high supersonic stage 
dvan der Swaluw et all 12004). This is consistent with an age of 
few thousand years, and in ro ugh agreement with th e mildly su- 
personic Mach value found by Gaensle r et al.l d2006l) . 

The orientation of the cometary tail and the position of the 
PWN cannot be straightforwardly explained. This mismatch is 
still an open issue, but it does not precludes the association be- 
tween the PWN, the ring and the SNR (see also lGaensler et all 
2006). For instance, a similar apparent offset between the stellar 
remnant, the forward, and the reverse shock geomet r ic cen ters 
has been already observed in Cas A ( Gotthel f et all 1200 ll) . A 
possibility is that strong reverse shocks, originated by the inter- 
action with the surrounding clouds, skewed the orientation of the 
nebula. An asymmetric SN explosion or the presence of strong 
density gradients across the remnant offer a valid explanatio n 
to the observed offsets dGotthelf et all 1200 UlReed et all 1 1995b . 
Another possible explanation is that the SN ejecta have been de- 
celerated by the southern ridge of molecular material, and, as 
a consequence, the geometric center of the ring appears shifted 
toward the north direction. 



5. Conclusions 

We analyzed the hard (1.4-5.0 keV) X-ray emission of the su- 
pernova remnant IC 443, as observed by a set of public archive 
XMM-Newton observations. Unlike the soft X-ray emission, 
which traces direct ly the interaction with the surrounding clouds 
dTroia et all [2006), the hard X-ray emission has a centrally- 
peaked morphology, reminiscent of the subclass of mixed mor- 
phology SNRs. The presence of a ring-shaped feature, which en- 
circles the PWN, has been revealed in the equivalent width maps 
of Si and S lines (Fig. 3). 

A spatially-resolved spectral analysis confirms that this 
structure is associated to hot metal rich plasma, whose abun- 
dances are consistent with a core-collapse SN origin, but not 
with an exploding white dwarf. We argue that the ring indicates 
the location of an ejecta layer heated by the reverse shock, and 
that the reverse shock has partially hit the PWN. Dynamically, 
the location of the ejecta ring suggests an age of ~4,000 yr for 
IC 443, an order of magnitude lower than some previous es- 
timates, based on expansion in a very dense environment. A 
younger age is indeed more consistent with an SNR expansion 
within wind-blown bub ble shells, surrounde d by dense clouds, 
as first su ggested bvlBraun & Stroml d!986l) . and recently con- 
sidered by Tr oia et all d2006l) . 



We detect only marginal evidence of plasma overionization 
in the NE part of the remnant. The bright hard 1.4-5.0 keV X- 
ray emission in the NE region of IC 443 also show high tempera - 
tures and high metal abundances, as shown in lTroia et al.l (|2006). 
More detailed results from this region will be reported in a future 
work. 
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